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Abstract 

We present a method for simulating the production and decay of particles in the sextet 
representation of SU(3)c including the simulation of QCD radiation. Results from 
the Monte Carlo simulation of sextet diquark production at the LHC including both 
resonant and pair production are presented. We include limits on resonant diquark 
production from recent ATLAS results and perform the first simulation studies of the 
less model dependent pair production mechanism. 



1 Introduction 

Many models of Beyond the Standard Model (BSM) physics require the inclusion 
of diquarks. For example, diquarks appear in a number of Grand Unified The- 
ories (GUTs) and have even been postulated as a form of dynamical symmetry 
breaking, giving rise to the masses of particles, [lj[2]- The colour sextet diquark 
is, in group theory language, a rank 2 symmetric tensor formed from the direct 
product of two fundamental representations 3 <S> 3 = 6 © 3. As such it is the 
lowest colour representation which has not been observed and therefore investi- 
gation of sextet diquark production at the CERN Large Hadron Collider (LHC) 
is interesting in its own right. 

The LHC experiments have started data taking at the high energy frontier 
(y/s = 7TeV), allowing probes of energy scales not previously seen. At the LHC 
because the fundamental collisions are between the quarks and gluons inside the 
colliding protons the strong force is the dominant interaction allowing Quantum 
Chromodynamics (QCD) to be studied at these new high energies. As a diquark 
is produced via strong interactions, and with the potential of a relatively low 
mass, diquarks may be seen in the early stages of LHC data taking. The LHC 
also favours the formation of diquarks from the valence quarks as it is a proton- 
proton collider as opposed to a proton-antiproton collider, such as the Tevatron. 

Due to their exotic colour structure and SU(3)c quantum numbers, diquarks 
will give rise to jets in the detector. The expected signals will be either a resonance 
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in the invariant dijet mass distribution or the production of two equal mass dijet 
systems in four jet events in the case of pair production. In order to study 
the experimental signatures of diquark production, a Monte Carlo simulation 
is required that includes the production of sextet particles, their perturbative 
decays and the full Monte Carlo machinery of showering (including the exotic 
colour structure) and hadronization. 

Although significant efforts have been made to study the resonant production 
of diquarks [3jj9] and also pair production [5 10 11 , a full study of experimental 



signatures including Monte Carlo simulations has not been performed. In this 
paper, we discuss: the implementation of the diquark model in general purpose 
Monte Carlo event generators; place constraints on the coupling as a function of 



mass based on the latest ATLAS results 12 13 ; present some results of invariant 



mass distributions both for resonant and pair production. 

2 Simulation 

Monte Carlo simulations describe high energy collisions using: 

1. a hard perturbative, either leading- or next-to-leading-order, matrix ele- 
ment to simulate the fundamental hard collision process; 

2. the parton shower algorithm which evolves from the scale of the hard process 
to a cut-off scale, 0(1 GeV), via the successive radiation of soft and collinear 
quarks and gluons; 

3. the generation of multiple perturbative scattering processes to simulate the 
underlying event; 

4. the perturbative decay of any fundamental particles, with lifetimes shorter 
than the timescale for hadron formation, followed again by the simulation of 
QCD radiation from the coloured decay products using the parton shower 
formalism; 

5. a multiple partonic collision model is used to simulate the underlying event; 

6. a hadronization model which describes the formation of hadrons at the cut- 
off scale from the quarks and gluons produced during the parton shower; 

7. the decays of the unstable hadrons produced by hadronization. 

The various calculations, approximations and models used in these simulations 



are reviewed in Ref. 14 . Simulating most models of BSM physics only requires 
the implementation of the various hard production and decay processes with the 
simulation of perturbative QCD radiation and hadronization proceeding in the 
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same way as in the Standard Model (SM), provided the new particles decay be- 
fore forming hadrons. However, in models including the production and decay 
of sextet particles or R-parity violating SUSY[j] the new colour structures require 
changes to both the simulation of QCD radiation and the subsequent hadroniza- 
tion. 

As with all models the simulation starts with the calculation of the hard pro- 
duction and decay processes using the most general Lagrangian for the coupling 
of the sextet particles to the quarks [4, 6,[9| [l6|[T7| 



+ (i) 



£ = {9iLq c L iT 2 qL + 9iRU c R d R ) $1,1/3 + g[ R d c R d R $i- 2 /z + 9irU r u r $ 1a / 3 
93Lqlir 2 Tq L ■ $3,1/3 + 92qJl^d R V£_ 1/6 + g 2 qJl,,u R V£ 5/6 + h.c. , 

where qi is the left-handed quark doublet, u R and d R are the right-handed quark 
singlet fields, and q c = Cq T is the charge conjugate quark field. The colour and 
generation indices are omitted to give a more compact notation and the subscripts 
on the scalar, $, and vector, V M , fields denote the SM electroweak gauge quantum 
numbers: (SU(2) L , U(l)y). The Lagrangian is assumed to be flavour diagonal 
to avoid any flavour changing currents arising from the new interactions. 
The kinetic and QCD terms in the Lagrangian are 

C s ^ r = D^Dfl - m 2 <$>&, (2a) 

for scalar diquarks, where $ is the scalar diquark field and 



£QCD r = ~\ (D^V U - D V V») (D,V U - D V V») - mV%, (2b) 



for vector diquarks, where is the vector diquark field. The covariant derivative 
has the standard form for Quantum Chromodynamics. 

The simulation of perturbative QCD radiation, relies on the large number of 
colours, Nc, limit for both the treatment of perturbative QCD radiation and the 
subsequent hadronization. In this approach particles in the fundamental repre- 
sentation of SU(Nc) carry a colour, those in the antifundamental representation 
an anticolour and those in the adjoint representation both a colour and an an- 
ticolour. This allows us to consider the flow of colour, via colour lines, in hard 
interactions which is determined by the colour structure of the hard perturbative 
matrix elements. 

Using this colour flow both: 

• soft gluon radiation, where in a given hard process there is a maximum 
angle for the radiation of gluon related to the colour flow of the process; 



and hadronization, where the first step of both the string 18, 19 and clus- 
ter [20] models is the formation of colour singlet systems; 



1 The simulation of R-parity violating SUSY models was considered in detail in Ref. 15 . 
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can be simulated. 

This is complicated in models involving sextet particles where in the large- 
Nq limit the sextet particles possess two fundamental colours, appropriately sym- 
metrized. This cannot be handled by conventional Monte Carlo simulations which 
require all the colours of the particles to have fundamental colours and/or anti- 
colours. In order to simulate these particles we choose to represent (anti)sextet 
particles as having two (anti) colours. 

Consider the production and subsequent decay of a scalar sextet particle. In 
order to simulate QCD radiation from the intermediate sextet resonance we have 
to simulate the production and decay separately. The matrix element for the 
process is 

M = M ipiod ^- 2 ML ay (3) 

where i,j are colour indices of the sextet particle, -Mj pro d is the matrix element 
for the production of a scalar sextet particle with colour i, four-momentum p and 
mass m, and A1^ ecay is the matrix element for the decay of a scalar sextet particle 
with colour j. 

This can be rewritten using £j = K l ab Kj a where K and K are the Clebsch- 
Gordan coefficients in the sextet and antisextet representations, respectively. 
Hence 

%K l K ba j 
M = M ^^^ M L, y = M' abpiod ^-^M>£ car (4) 

In order to consider the intermediate sextet particle as having two fundamen- 
tal colours we have absorbed the Clebsch-Gordan into the redefined production, 
■M.' a bp ro d, and decay matrix elements, -M' ^ cay . 

From this, the colour partners of the decay products can be determined and 



the usual angular ordering procedure applied 14,21 . The radiation pattern of 
gluons from the quarks for the resonant production of diquarks, 
q(Pi)q(P2) -> {^,V} q(p 3 )q(Pi), is 



J 2 = C f 



{W 13 + W 14 + W 23 + W 2A ) (5) 

\{w l3 + w u + w 23 + W 24 ) - w 12 - w 34 



+ 



N c + l 

The dipole radiation function for two massless particles i and j radiating a gluon 
is 

where Pij are the 4-momenta of the radiating particles and, k and u are the 
4-momentum of energy of the gluon, respectively. In terms of the angle between 
the gluon and particle i, 9i, and the angle between the particles i and j, dij, 



W*- = 1 fl + C ° S ^~ CQS% ) (7) 

l] 2(1- cos OA V 1- cos ft ' ' 
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Fig. 1: The radiation pattern of associated with gluon emission from the incoming 
and outgoing quarks during resonant production, where 9 is the polar 
angle of the gluon with respect to the z-axis. 



The last term in Eq. 5| can be neglected, as usual, due to both the sup- 
pression, compared to the leading term, and the dynamical suppression in the 
massless limit because there is no collinear singularity in this term. The radiation 
pattern is shown in Figure [TJ where the outgoing quarks were held at 45° and 
225° with respect to the incoming beam direction. The full radiation pattern, the 
result after neglecting the subleading terms and azimuthally averaging, and the 
improved angular ordered result, where the full result is used instead of the az- 
imuthal average inside the angular-ordered region, are shown. Improved angular 
ordering, as implemented in Herwig-r-+ performs well in the collinear limit. 

The details of the colour decomposition for both the resonant and pair pro- 
duction of sextet particles and the radiation of gluons from sextet particles are 
presented in Appendix A. 3 The same approach was recently used in Mad- 
Graph 5 22 including the ability to automatically decompose the colours for 



higher representations of the gauge group. They also study the impact of match- 
ing the hardest perturbative emission in resonance production which we have not 
considered but will not effect on our results. 



3 Phenomenology 

In order to study the phenomenology simulations were performed for the scalar 
$1,4/3 and the vector V^ S / 6 diquarks which were chosen as they can be produced 
as s-channel resonances from the partonic collision of the valence up quarks. 

, as 



In all our analyses, jets were clustered using the anti-fcr algorithm 23 
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Mass [GeV] 

(a) Scalar 




Mass [GeV] 

(b) Vector 



Fig. 2: The Breit-Wigner shapes for y/s 
and 0.05. 



14TeV with couplings of 0.1, 0.075 




Mass [GeV] 



Fig. 3: The diquark width as a function of the diquark mass for a coupling as 
quoted in the text for scalar and vector diquarks. The diquark coupling 
to quarks has been taken to be 0.1. 



implemented in the FastJet package 24 , using a radius parameter of R = 0.6. 
This choice is typical for the ATLAS experiment at the LHC 12 . The LO** 



PDFs of Ref. [25], which are the default choice in Herwig++, were used. 

There are phenomenological constraints on the diquark couplings from D° — D° 
mixing and non-strange pion decays 17]. For the up-type quarks there are con- 
straints require that 

g u R u <0.1 and g%~0. ( 8 ) 
The gi couplings have to be constrained due to minimal flavour violation as the 
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Fig. 4: Resonant production and decay of a diquark 



left-handed CKM matrix is well known (31. 

It was therefore decided to take the couplings 

9$/l = M an d 9%l = 9%l = ^ (9) 

where the numbered indices refer to the generation. 

The value of the coupling will affect any studies involving jets as the width 
of the particle varies as a function of the couplings. As seen in Figure |2j a larger 
coupling produces a larger width, contributing to the smearing of a peak in the 
invariant mass spectrum of the jets. If the coupling is less than the value chosen 
above, any peak maybe enhanced compared to what is presented in the following 
sections. The width as a function of the diquark mass is shown in Figure |3j 



3.1 Resonance Production 

If the diquark has an appropriate mass and coupling it may be resonantly pro- 
duced at the LHC. The resonance production and subsequent decay of a general 
diquark (scalar or vector) is shown in Figure |4j where the decay of the diquark 
will depend on its mass and unknown couplings to the quarks, g. 

Figure [5] shows the cross section for the production of scalar and vector di- 
quarks from incoming uu quarks (resonant production) and for incoming gluons 
(pair production). The resonant production cross section depends quadratically 
on the unknown diquark coupling to quarks, which as been assumed to be 0.1 
in this plot, whereas the pair production cross section is independent of this 
coupling. 

The diquark will decay into two quarks giving rise to at least two jets. The 
search for the production of a diquark via resonant production should therefore 
be in the dijet invariant mass spectrum, where a smeared peak is expected around 
the diquark mass. The primary background to this search channel is QCD 2—7-2 
scattering processes. 

The signal and background were simulated using Herwig+-|-. The analysis 



and modelling of the backgrounds followed that suggested in Ref. 12 . The 



transverse momenta and pseudorapidities of the jets were required to be p\ > 
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Fig. 5: Cross section for the production of vector and scalar diquarks as a func- 
tion of the diquark mass for both resonant production, from incoming uu 
states, and diquark pair production at y/s = 14 TeV. The diquark coupling 
to quarks has been taken to be 0.1. 



150 GeV, p\ > 60GeV and |^i )2 | < 2.5 where 1 is the hardest jet and 2 is the 
sub leading jet. In addition we required that the dijet invariant mass, rrijj satisfied 
rrijj > 300 GeV and the rapidity difference between the leading and subleading 
jet was IA7712I < 1.3. The dijet invariant mass spectrum after these cuts is shown 
in Figures [6] and [7] for y/s = 7 and 14 TeV, respectively. The diquarks were 
simulated at masses of 500 GeV, 800 GeV, 1200 GeV, 1600 GeV and 2000 GeV. 

As simulating the QCD rrijj spectrum at high masses is difficult, a functional 
form 

f(x) = a Q (l - x)x {ai+a2lnx \ (10) 

was fitted to the low masses, where the are fitted parameters and x = rrijj /y/s, 
and extrapolated out into the high mass region. 



The results of Ref. 12 can be used to impose constraints on the diquark 

was 



coupling as a function of the diquark mass. The event selection from Ref. 12 
used to reproduce the correct acceptance. This requires that the event contains 
at least two jets with px > 150 GeV and a subleading jet with px > 30 GeV. 
Both the leading and subleading p T jets must satisfy \r]j\ < 2.5 with Ar] 12 < 1.3 
and rrijj > 150 GeV. 

The signal, after the cuts, was fitted to a Gaussian distribution, with the mean 
fixed, m, at the simulated diquark mass to obtain the standard deviation of the 



distribution, a, so that the results presented in Ref. 12 could be used to obtain 
the limits on the diquark coupling. As suggested in |12| , long tails were removed 
by taking a window around the diquark mass of ±20% for the fit. If the a/m 
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Reconstructed Mass mjj [GeV] Reconstructed Mass mjj [GeV] 



(a) Scalar (b) Vector 



Fig. 7: The dijet mass spectrum at ^ = 14 TeV for 500 GeV, 800 GeV, 1200 GeV, 
1600 GeV and 2000 GeV diquark masses with the couplings given in the 
text. 

value obtained was below the range of that given in the paper, then the number 
of events associated with the lowest a/m for that mass was used. This allows 
a conservative estimate for the excluded coupling, as opposed to one which may 
be obtained by extrapolation into the unknown region. The limit on the diquark 
coupling is shown in Figure [8] where because the statistical errors were negligible, 
the bands shown come from varying the scale from 50% to 200% of the default 
scale choice, i.e. the diquark mass. 



Following the work of Ref. 12 the ATLAS collaboration has released an 
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£ = 36pb- 1 

- £ = 163pb- 1 
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Fig. 8: Limit on the coupling as a function of the diquark mass based on the 



model independent data given in the recent ATLAS publications 12 , 13 
The band reflects the uncertainty from varying the scale between 50% and 
200% of the diquark mass. 



updated analysis 13 , including additional data corresponding to an integrated 
luminosity of 163 pb -1 . This analysis included slightly harder cuts requiring 
Pt > 180 GeV and rrijj > 170 GeV in addition to the cuts used in Ref. 12 



The limit obtained from this higher integrated luminosity analysis is also shown 
in Figure |8l We note that ATLAS performed better than the expected median 



limit in the Ref. 12 and worse than the expected median limit in 13 in the 



1400 — 1600 GeV mass range, giving rise to the overlap in Figure [8j 

3.2 Pair Production 

The pair production of diquarks (scalar and vector) occurs via the Feynman di- 
agrams shown in Figure [9] The pair production process has one main advantage 
over the resonant production, it does not depend on the unknown diquark cou- 
pling. Instead, the pair production cross section depends only on the SU(3)c 
representation, mass and spin of the particle. The pair production processes 
therefore has the potential to distinguish between whether a particle in the an- 
titriplet or sextet representation was produced due to the dependency of the cross 
section on the colour representation. 

To date there have been no studies of the experimental signals of diquark 
pair production. The cross section has been calculated [5, 10 11 and some work 



towards a jet study, no Monte Carlo study has been performed. 

The pair production and subsequent decay of diquarks is expected to give four 
jets, with two pairs of jets forming systems with the mass of the diquark. The 
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(a) s-channel 






(b) f-channcl 




(c) u-channel (d) 4-point Channel 

Fig. 9: The four diagrams contributing to the pair production of a diquark. 



backgrounds to the pair production of diquarks are: 

• vector boson WW, ZZ and ZW pair production; 

• vector boson, W and Z, production in association with additional jets; 

• top quark pair, ti, production; 

• QCD jet production. 

The analysis proceeded by placing cuts on the four hardest jets: p\ > 150 GeV, 
■p\ > 100 GeV, > 60 GeV and p\ > 30 GeV, where the four jets i = 1,4 are 
ordered in p? such that the first jet is the hardest. All four jets were required to 
have pseudorapidity \r)i\ < 3. Two pairs of jets were then formed with the pairing 
selected that minimized the mass difference between the two pairs of jets. If, after 
pairing, the two hardest jets are in the same pair of jets, the event was vetoed. 
The mass difference between the pairs was required to be less than 20 GeV. 

The signal and backgrounds were simulated for the production of the $1,1/3 



10 



and 



E 



for 



and V^ 5 , 6 diquarks giving the results shown in Figures 
y/s = 7 and 14TeV, respectively. As the backgrounds are dominated by QCD 
scattering, i.e. the contribution of the QCD scattering processes is approxi- 
mately one hundred times that of all the other backgrounds combined, only the 
sum of the backgrounds is shown. As for the resonant production, the plots 
show the production of diquarks with masses of 500 GeV, 800 GeV and 1200 GeV, 
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Fig. 10: The mass spectrum of dijet pairs in four jet events at y/s = 7TeV for 
500 GeV, 800 GeV, 1200 GeV, 1600 GeV and 2000 GeV diquark masses. 




Fig. 11: The mass spectrum of dijet pairs in four jet events at y/s = 14TeV for 
500 GeV, 800 GeV, 1200 GeV, 1600 GeV and 2000 GeV diquark masses. 



1600 GeV and 2000 GeV. The low mass QCD background was fitted with Eq. [10 
and extended out into the high mass region. 

A window of ±50 GeV was taken around the diquark mass and the ^= was 



calculated for a number of luminosities the results of which are shown in Figure 12 



We see that that the vector diquark manifests itself more prominently than the 
scalar, which is consistent with the increased cross section of the vector in the 
pair production process as seen in Figure |5j It will be hard to observe a scalar 
diquark using the pair production process at y/s = 7TeV while with £ = 5fb _1 
it should be possible to observe a vector diquark with mass less than 700 GeV. 
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Fig. 12: -7= for the scalar and vector diquark at luminosities of £ = lfb 1 
and C = 5ft- 1 at = 7TeV and C = lOfb" 1 and C = lOOfb" 1 at 
y/s = 14TeV. The black horizontal line shows = 5. 

There is a marked increase in discovery potential at the increased energy and 
luminosities running at y/s = 14 TeV brings. A vector diquark in the mass range 
presented here (< 2000 GeV) should be seen with C = 10 fb -1 , whereas even 
with C = 100 ftT 1 only a low mass (< 1050 GeV) scalar diquark has potential for 
discovery 

4 Conclusions 

In this paper we have presented a method for simulating the production and 
decay of particles in the sextet colour representation. This approach has been 
implemented in Herwig++ and will be available in a forthcoming release. 

This new simulation was used to simulate the production and decay of vector 
and scalar sextet diquarks at energies relevant to the LHC Based on the find- 
ings and the latest ATLAS search for new particles in two-jet final states, new 
constraints have been put on the couplings of the diquarks to SM particles. 

We have presented the first studies of the pair production mechanism, which 
is independent of the unknown coupling of the sextet diquark to quarks. This 
process has a promising search reach with the possibility of observing vector 
diquarks with masses less than 710 GeV at y/s = 7 TeV and both vector, for 
masses less than 2 TeV, and scalar, for masses less than 1 TeV, with the LHC 
running at design energy. Hopefully the availability of a Monte Carlo simulation 
of these processes will allow a more detailed experimental study. 
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A Colour Decomposition 

Herwig++, as with all general purpose Monte Carlo generators, has all the ma- 
chinery set up to work with the (anti)fundamental representation of SU(3)c- 
Therefore exotic color representations must be decomposed into a fundamental 
representation basis. 

The results of the this decomposition for the resonant and pair production of 
sextet diquarks and the radiation of a gluon from a sextet diquark are explicitly 
outlined below. 

A.l Resonant Production 

The diagram for resonant production and decay is shown in Figure |4} Following 
decomposition into the fundamental representation, there are two unique colour 



flows associated with this process as shown in Figure 13. The colour factor 



N C (N C +1) 



associated with these colour flows is 





(a) (b) 
Fig. 13: Unique colour flows associated with the resonant production of a diquark. 



A. 2 Pair Production 

The diagrams contributing to the pair production of diquarks are shown in Fig- 
ure |9j Each of these diagrams is decomposed into a colour factor and colourless 
component. By taking the colour factors for the diagrams, adding and squaring, a 
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table of colour factors can be produced for each term in the total matrix element 
squared. 

The pair production process has twelve unique colour flows, as shown in Fig- 
ure 14 The colours in Figure 14 have no physical meaning and are included as a 



visual aid. 

From these colour flows, Table 1 is produced, where c\ 

c 3 = ^f^, c 4 = ^fe?, c 5 = ^f 11 and c 6 = 



-i) 2 



16N C 



16N C 



16 



1) 



16N C 



-l) 2 



16 
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Tab. 1: Associated colour factors for the diagrams shown in Figure 14 The values 
of q are given in the text. 



A. 3 Shower 

The splitting functions were decomposed in the same way as for the pair pro- 
duction. There are four unique colour flows associated with a diquark emitting 



a gluon, as shown in Figure [15j again where the colours are included as a visual 
aid. 

Only the colour prefactor of the existing splitting functions is changed. The 
colour prefactor is given by y, i.e. the diquarks radiate 2| times more than a 
particle in the octet representation. 

During the shower, it is assumed that gluons did not branch to form diquarks 
owing to the large diquark mass. 
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0) (k) (1) 

Fig. 14: Unique colour flows associated with the pair production of diquarks. 
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(a) (b) 



(c) (d) 
Fig. 15: Colour flows for a diquark emitting a gluon during the shower. 
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